Polar decomposition and Brion's theorem. 



Christian Haase 

Abstract. In this note we point out the relation between Brion's 
formula for the lattice point generating function of a convex poly- 
tope in terms of the vertex cones [Bri88] on the one hand, and the 
polar decomposition a la Lawrence/ Varchcnko [Law91a, Var87] 
on the other. We then go on to prove a version of polar decompo- 
sition for non-simple polytopes. 



1. Introduction 

The stars of this note are two formulas that express convex polytopes 
in terms of cones. On the one hand, Brion's theorem (equation (2) 
below) expresses lattice points in a polytope in terms of tangent cones 
at vertices [Bri88]. On the other hand, given a simple polytope and 
a generic objective function, one can write the indicator function of 
a simple polytope as a combination of polarized tangent cones at the 
vertices (equation (4)). We refer to this formula as polar decomposi- 
tion [Law91a, Var87]. 

After we review these formulas in this introduction, we point out in 
section 2 some simple relations between the two formulas which, to our 
knowledge, have not appeared in the literature before. In particular, 
Brion's theorem (for simple polytopes) follows from polar decomposi- 
tion. The last section is devoted to the formulation and proof of polar 
decomposition for non-simple polytopes 1 . 

This is by no means an attempt to survey what is scattered through- 
out the literature 2 about decompositions of polytopes into cones. In 

Work supported by NSF grant DMS-0200740. 

x Put together, sections 2 and 3 yield another proof of Brion's theorem. 

2 See [Bar02, p. 346] and the references therein for Brion's theorem, [Law91a, 
Var87] for polar decomposition, and [Bri37, Gra74, She67] for the Brianchon- 
Gram formula below. 
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fact, Matthias Beck made me aware of such a survey in the mak- 
ing [BRS W04] . We will concentrate on what we think is new. We will 
not state the most general form of the results, but leave generalizations 
(e.g., weighted versions) as exercises for the ambitious reader. 

Acknowledgments. This paper has its origins in Jose Agapito's in- 
spiring talk about his weighted version of polar decomposition for sim- 
ple polytopes [Aga03]. After the talk, I asked the question whether 
one can obtain an analogous formula for non-simple polytopes by con- 
sidering simple deformations. I got intrigued by this question, and, 
after computing several examples, was convinced that there is such a 
theorem out there. 

I was lucky to get hold of a preliminary draft of the above mentioned 
survey [BRSW04]. Not only did it take away the pressure to try and 
write a comprehensive article about conic decompositions of polytopes, 
but it also made me aware of references and view points that I did not 
know about before. Much of the presentation in the present paper is 
influenced by (if not stolen from) this draft. Finally, I want to thank 
the referee who helped to clarify the exposition. 

1.1. Basic definitions and the Brianchon-Gram formula. 

We will consider rational convex polytopes P and polyhedral cones C 
in Q d , where Q d is endowed with the fixed lattice Z d . We will assume 
that polytopes and cones are full-dimensional. For standard polytope 
definitions and notation we refer to [Zie95]. To a subset S C Q d we 
assign two objects: 

- The indicator function t s : Q d — > Z vanishes outside S, and is one 
along S. 

- The generating function (of the lattice points) is the formal Laurent 

power series G s = J2mesnz d e ^H 2 ! 11 ' • • • > z d *]]> where we write 
z m = z™ 1 ■ . . . ■ z™ d . If S = P is a polytope, then Gp is a Laurent 
polynomial; when S = C is a cone, then Gq is a rational function 3 . 

The essential features can be illustrated on the line. For a < b G Z, we 
have G [aM (x) = x a + . . . + x b = xa ~f^ , an d G> )0o) (x) = x a + ... = 

The mother of all conic decomposition theorems is the Brianchon- 
Gram formula. According to Shephard [She67], Brianchon [Bri37] 
and Gram [Gra74] independently proved the d — 3 case in 1837(!) and 
1874 respectively. In 1927 Sommerville [Som27] published a proof for 

3 We are glossing over significant technicalities here and later on when we pass 
from formal power series to rational functions. For a first reading, we prefer to skip 
this step, and refer to [Bar02, BV97] for a careful treatment. The experts should 
be familiar with the standard arguments, anyway. 
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general d, which was corrected in the 1960's by Griinbaum [Grii03, 
§14.1]. 

Let P C Q d be a convex polytope. For a face F of P, we define the 
tangent cone (of P at F) by 

T F P = {/ + x G Q d : / G F and / + ex G P for some £ > 0} 

(See Figure 1.) Then lp is the alternating sum of all tangent cones: 




Figure 1. Tangent cones at vertex v, edge e, and at P. 



(1) lp=^(-l) dim ^lr F P. 

F^P 

At every point x G F, the right hand side computes the Euler char- 
acteristic of F (and F is contractible). While outside F, we have to 
subtract the Euler characteristic of the subcomplex that is visible from 
x (which again, is contractible, actually shellable). (Compare [She67].) 

1.2. Brion's formula. Brion [Bri88] proved his decomposition 
using a Riemann-Roch type formula. In the mean time, more elemen- 
tary proofs have been given [Bar93, BV97, Law91b, PK92]. The 

easiest way to formulate Brion's formula is in terms of generating func- 
tions. We can express the generating function for F as the sum of the 
generating functions of the tangent cones at vertices. 

(2) G P = G ^ 

v vertex of P 

In the example, T a [a,b] = [a, oo), and %[a, b] = (—00,6] with gen- 
erating functions G [cit0o) = ^ and G(_oo,&] = = -f=j- As 
predicted, they sum to G[ a ,&]- 

This version of Brion's formula is implied by the following two re- 
sults. (Compare [Bar02].) 

Theorem 1. The function t P — vertex j P t%p is a linear combi- 
nation of indicator functions of cones that contain affine lines. 

This follows immediately from (1). It is also implied by the consid- 
erations in Sections 2 and 3. In our example, l[ a ,b] — l[a,oo) + l(oo,b] — 1q- 
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LEMMA 2. Let C C Q d be a cone that contains an affine line. Then 
G c = 0- 

The idea of the proof is to decompose the monoid C H Z, d into a 
direct sum L © R of a one dimensional lattice L and a complement R. 
Then Gc can be written as GrGl, and Gl = 4 . 

1.3. (Simple) polar decomposition. By flipping the edge vec- 
tors emanating from each vertex of P in a systematic way, we can get 
the so-called polar decomposition, which expresses the characteristic 
function of a convex simple polytope in terms of the characteristic func- 
tions of the cones supported at the vertices of P (no cones with straight 
lines needed this time). This was first obtained by Varchenko [Var87] 
and Lawrence [Law91a] independently. Lawrence used Brianchon- 
Gram's theorem together with the principle of inclusion-exclusion to 
derive the polar decomposition theorem. Then Karshon, Sternberg 
and Weitsman obtained a weighted version of this decomposition by 
assigning, in a consistent way, particular weights to the lattice points 
in the polytope and in the cones [KSW03]. Their work was also moti- 
vated by methods in differential and algebraic geometry. Finally, using 
the same source of motivation, Agapito [Aga03] gave a more gen- 
eral weighted version of the polar decomposition theorem for simple 
polytopes that includes Lawrence/ Varchenko and Karshon-Sternberg- 
Weitsman versions as particular cases. 

Polar decomposition is like Morse theory for simple polytopes. We 
sweep over the polytope and build it up from local contributions, crit- 
ical point (vertex) by critical point. As input we use a generic 5 linear 
functional £ G (Q rf )*. Using £, we define the polarized tangent cones 
as follows. 

Because P is simple, at every vertex v, the tangent cone T V P is 
generated by d linearly independent directions t±, . . . ,td (remember P 
is full-dimensional). 

d 

%P = v + J2^>o^ 

i=l 

Now the polarized tangent cone (with respect to £) is 

Q<oU- 

«(*i)>0 ?(ti)<0 

This is a locally closed cone all whose points are ^-higher than v. 
4 Glossing over alert: look under the rug [Bar02]! 

5 Generic means here that £ is not constant on any edge of P. In Section 3 we 
ask a little bit more. 
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FIGURE 2. Polarized tangent cones. 

(Compare Figure 2.) Dually, if n\, . . . ,nd G (Q rf )* are the inner facet 
normals to P at t), and £ = c^rii, then 7;«P can be defined by the 
inequalities 

(3) rii(x) > riiiv) for on > 0, and rii(x) < rii(v) for on < 0. 

Define the index ind^(u) of the vertex v as the number of ^-negative 
edge directions ij, or equivalently, as the number of negative coefficients 
«j. Then we can write P as the signed sum of polarized tangent cones. 

(4) i P = Yl (-i) in<w %p 

v vertex of P 

Again, we can formulate a generating function version of this result: 
just replace all 1's by G"s 6 . 

Equation (4) follows from (1) if we group all faces according to 
where they achieve their ^-maximum. Then it remains to check that 

(5) (-ir d ^i %(p = Yl (-i) dimF it>p. 

v<F<P 

There is a weighted version of equation (4) [Aga03]. The weighted 
indicator function of a polytope P is the function If, : Q d — > Z[z] which 
vanishes outside P, and takes the value z k along the relative interior 
of codimension k faces. 7 One could even introduce a different variable 
for every facet of P, and assign their product to the corresponding 
intersection. Now the same formula holds if we modify the indicator 
functions of the polarized tangent cones as follows: A face of (the 
closure of) T^P is defined as the set of points that satisfy equality for 
some of the inequalities (3). The value of 1?S along a face defined 

by k + equalities rii(x) = rii{v) for «j > and A;_ equalities for ccj < 
is the polynomial z k+ {l — z) k - G Z[z]. One recovers the unweighted 
version for z = 1, and for z = one obtains a decomposition of the 
interior of P. 



'More glossing over. 

The formula in [Aga03] is stated using the substitution z = j^- . 
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2. This follows from that 

In this section we show the relation between the two formulas. Essen- 
tially, for the indicator functions, Brion's formula (in the simple case) 
is polar decomposition modulo cones that contain lines. 

Lemma 3. If P is a simple polytope and v is a vertex, then the differ- 
ent polarized tangent cones for various £ partition Q d . Moreover, for 
6,6 e (Q d Y, (-l) hld ^ {v) l^ lp and (-l)H a Wi, are equivalent 
modulo cones that contain lines. 

Proof. If rij(x) > riiiv) for % = 1, . . . , d are the inequalities of the 
facets incident to v, the hyperplanes Ui(x) = Ui(v) subdivide Q d into 
orthant cones. Every point x G Q d belongs to exactly one polarized 
tangent cone according to the signs of the nj(x). 

We can get from every polarized tangent cone to every other by suc- 
cessively flipping inequalities. Now the sum of two adjacent polarized 
tangent cones is defined by d — 1 (strict and non-strict) inequalities. 
Thus it is a cone that contains a line. □ 

For a simple polytope, the indicator function version of Brion's 
formula (Theorem 1) can be derived from equation (4) (modulo cones 
that contain lines). As Theorem 1 does not specify these cones, there 
is no converse. On the other hand, generating functions do not see 
cones that contain lines. So for generating functions, the two formulas 
actually coincide. 

The same considerations carry through for the weighted indicator 
functions. One thus obtains a weighted Brion's formula. We leave this 
as an exercise for the reader. 

3. Non-simple polar decomposition 

We want to generalize polar decomposition to non-simple polytopes. 
We will compute the local contribution from a simple deformation of 
the vertex in question, and prove that the result does not depend on 
the chosen deformation. 

As there are not too many non-simple polytopes in dimension < 2, 
the new running example will be the pyramid with the five vertices 
(0,0,0), and (±1,±1,1). 

Let ni, . . . , riN be the inner facet normals to P at u. They generate 
the rays of the normal cone M v . A virtual deformation of the vertex v is 
a regular triangulation 8 A of Af v . That is a face-to-face subdivision of 

8 Compare [Lee97, § 14.3]. Moving the facets according to the values of the 
convex function at the n.; would yield an actual simple deformation of the vertex v. 
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Figure 3. The simplest non-simple polytope. 

M v into simplicial cones a±, . . . , o~m , so that there is a convex piecewise 
linear function on M v with domains of linearity the Oi. 




Figure 4. Triangulations and the corresponding deformations. 

In the example, v = (0, 0, 0) is the only non-simple vertex. The nor- 
mal cone is generated by ri\ = (1, 0, 1), n 2 = (—1, 0, 1), n 3 = (0, 1,1), 
and ri4 = (0, —1, 1). It has two triangulations, both of which are reg- 
ular: Ai with maximal simplices spanned by 77,1,713,714 and 712,713,714, 
and A2 with maximal simplices spanned by 711,77.2,713 and 711,712,714. 
The corresponding deformations are sketched in Figure 4. 

Now the tangent cone T V P can be written as the intersection of 
simple cones defined by the a^. 

(6) %P = f]% { for T a . = {x EQ d : n^x) > n^v) for all n,j E a { }. 

i 

For generic 9 £ G (Q d )*, we write £ in terms of the 77j G CTj, and flip those 
inequalities in (6) where £ has negative coefficients. Then we compute 
the index ind^Uj) as the number of flipped inequalities, and define the 
local contribution at v as: 

lL = E(- 1 ) md5( " ) %r CT! 

i 

The cool thing is that 1^ v = l£ does not depend on the triangulation. 
So it really is a local contribution. 

9 Now, generic means that £ is not constant on any ray of any of the 7^'s, i.e., 
£ does not lie on any hyperplane used in the triangulation. 
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Theorem 4. Let P c Q d be any poly tope, and let£ G (Q d )* be generic. 
Then 1 ^ „ = if does not depend on the choice of a regular triangulation 
A of the inner normal cone J\f v at vertex v. Moreover, 

(7) l P (x)= Yl ^ 

v vertex of P 

The A-invariance can be shown, e.g., using the fact that all regu- 
lar triangulations are connected by flips [Lee97, § 14.6]. In this note, 
however, we follow a different strategy. First, we will show the decom- 
position formula (7) for compatible choices of triangulations (Lemma 
5). Then, we observe that under relatively weak conditions, all decom- 
positions are the same (Lemma 6). 

But let us first see how the A-invariance works out in our example. 
If we choose £ = (4,2,0), then we compute for triangulation Ai as 
follows: £ = 4ni — n 3 — 3n 4 = — 4n 2 + 3n 3 + n 4 . Thus the polarized 
cones are 

T^T ai = {x : ni(x) > 0,n 3 (x) < 0,n 4 (x) < 0} with index 2 and 

T,^T a2 = {x : n 2 (x) < 0,n 3 (x) > 0,n 4 (x) > 0} with index 1. 

The trace of 1^^- — I^t- in the £ = 12 plane is sketched in Figure 5 
on the left. Similarly, for A2 we get 



(6,-6,6) (6,-6,6) 




(6,-6,-6) (2,2,-2) (6,-6,-6) (2,2,-2) (0,6,0) 

Figure 5. The local contribution from Ai and from A2. 

TyT(T3 — { x '■ TT'x(x) > 0,ri2(x) < 0,n 3 (x) > 0} with index 1 and 

T^T a4 = {x : nx(x) > 0,n 2 (x) < 0,n 4 (x) < 0} with index 2. 

The trace of 1^^- — l^o- in the £ = 12 plane is sketched in Figure 5 

on the right. Observe that T}T az and T}T ai cancel out on the overlap 
to yield the same contribution as 1^T ai and T^T a2 . 
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Lemma 5. Let P C Q d be any polytope, and let £ G (Q d )* be generic. 
Suppose A is a regular triangulation of the polar dual polytope P* . It 
restricts to regular triangulations A v of the normal cones M v . Then 

(8) lp(x)= ^o*) 

v vertex of P 

One can mimic the proof for the simple version (4) in order to show 
that the sum of the 1^ v over all vertices of P does not depend on £, 
and that for every x G Q d there is a suitable £ such that equation (8) 
is satisfied. 

Now we actually use (8) together with the following lemma in order 
to show that 1^ v — l| is independent of the triangulation. 

Lemma 6. Let P c Q d be any polytope, and let £ be a generic element 
o/(Q )*. Suppose that there are two decompositions 

ip = Y f*> = S 9v 

v vertex of P v vertex of P 

that are both conic: for every direction t G Q d \ 0, the f v (v + \t) and 
g v (v + Xt) are constant in A > 0, and positive: f v (v + t) — g v (v + t) — 
if€(i) < 0. Then f v = g v . 




Figure 6. Positive conic decomposition of Ip is unique. 

Proof. This proof illustrates what was meant by Morse theory 
earlier on. We sweep over the polytope bottom to top, and only at the 
vertices does something happen. That is, we compare the restrictions 
of the functions to the hyperplanes £ = c. Order the vertices v±, . . . , Vk 
in ^-increasing order. For c G £(ui+i)[, we have 

k i—1 i—1 k 

^ = Yl 9v 3 = i P -j2 9v 3 = ip-J2 A = Yl - ^ 

j=i j=l j=l j=i 

along £ = c. For (p) we use positivity of the f v 's and ^„'s, and for (c) 
we use induction, and the fact that the /„'s and g v 's are conic. □ 
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Another incarnation of the positive conic decomposition is to group 
the summands in the Brianchon-Gram formula (1) as we did in equa- 
tion (5), though it is less obvious that this is positive. 

3.1. Homework, or Why this is no good. This is the wild 
speculation section. At the risk of exposing the full extent of my igno- 
rance, I ask a bunch of questions that I stumbled over while compiling 
these notes. 

If a good proof is one that makes us wiser [Man77], then this is 
a bad proof. It would be nicer to have a deformation independent 
definition of the local contribution in terms of the facet hyperplane 
arrangement at each vertex. How much geometry is needed? Can we 
compute the local contribution in purely combinatorial terms? 

In some sense, the actual deformation of a vertex into several simple 
vertices converges to the non-simple vertex, and the local contribution 
is continuous. While this asks for a topology on the polytope algebra, 
I believe that this is a more combinatorial question. There must be a 
combinatorial deformation theory for (Eulerian?) posets. 

Where are Grobner bases? Often when regular triangulations show 
up, the corresponding convex function determines a term order. Does 
the result tell us something about commutative algebra? 



References 

[Aga03] Jose Agapito. A weighted version of quantization commutes with re- 
duction principle for a toric manifold, July 2003. Preprint math.SG/ 
0307318. 

[Bar93] Alexander I. Barvinok. Computing the volume, counting integral points, 
and exponential sums. Discrete Comput. Geom., 10(2):123-141, 1993. 

[Bar02] Alexander Barvinok. A course in convexity, volume 54 of Graduate 
Studies in Mathematics. American Mathematical Society, Providence, 
PJ, 2002. 

[Bri37] C. J. Brianchon. Theoreme nouveau sur les polyedres. J. Ecole Poly- 
technique, 15:317-319, 1837. 

[Bri88] Michel Brion. Points entiers dans les polyedres convexes. Ann. Sci. Ecole 
Norm. Sup. (4), 21(4):653-663, 1988. 

[BRSW04] Matthias Beck, Sinai Robins, Frank Sottile, and Jonathan Weitsman. 
Conic decompositions of polytopcs, 2004. In preparation. 

[BV97] Michel Brion and Michele Vcrgne. Residue formulae, vector partition 
functions and lattice points in rational polytopcs. J. Amer. Math. Soc, 
10(4): 797-833, 1997. 

[Gra74] J. P. Gram. Om rumvinklerne i et polyeder. Tidsskrift for Math. ( Copen- 
hagen), 4(3):161-163, 1874. 



POLAR DECOMPOSITION AND BRION'S THEOREM. 



11 



[Gru03] Branko Grunbaum. Convex polytopes, volume 221 of Graduate Texts 
in Mathematics. Springer- Verlag, New York, second edition, 2003. Pre- 
pared and with a preface by Volker Kaibel, Victor Klec and Giinter M. 
Ziegler. 

[KSW03] Yael Karshon, Shlomo Sternberg, and Jonathan Wcitsman. The Eulcr- 
Maclaurin formula for simple integral polytopes. Proc. Natl. Acad. Sci. 
USA, 100(2):426-433 (electronic), 2003. 

[Law91a] Jim Lawrence. Polytope volume computation. Math. Comp., 
57(195):259-271, 1991. 

[Law91b] Jim Lawrence. Rational-function-valued valuations on polyhedra. 

In Discrete and computational geometry (New Brunswick, NJ, 
1989/1990), volume 6 of DIMACS Ser. Discrete Math. Theoret. Corn- 
put. Sci., pages 199-208. Amer. Math. Soc, Providence, RI, 1991. 

[Lee97] Carl W. Lee. Subdivisions and triangulations of polytopes. In Jacob E. 

Goodman and Joseph O'Rourke, editors, Handbook of Discrete and 
Computational Geometry, pages 271-290. CRC-Press, New York, 1997. 

[Man77] Yu. I. Manin. A course in mathematical logic. Springer- Verlag, New 
York, 1977. Translated from the Russian by Neal Koblitz, Graduate 
Texts in Mathematics, Vol. 53. 

[PK92] A. V. Pukhlikov and A. G. Khovanskh. The Riemann-Roch theorem for 
integrals and sums of quasipolynomials on virtual polytopes. Algebra i 
Analiz, 4(4):188-216, 1992. 

[She67] G. C. Shephard. An elementary proof of Gram's theorem for convex 
polytopes. Canad. J. Math., 19:1214-1217, 1967. 

[Som27] D. M. Y. Sommerville. The relations connecting the angle-sums and 
volume of a polytope in space of n dimensions. Proceedings Royal Soc. 
London (A), 115:103-119, 1927. 

[Var87] A. N. Varchenko. Combinatorics and topology of the arrangement of 
affine hypcrplanes in the real space. Funktsional. Anal, i Prilozhen., 
21(l):ll-22, 1987. 

[Zie95] Giinter M. Ziegler. Lectures on Polytopes, volume 152 of GTM. 

Springer- Verlag, 1995. 
E-mail address: haase@math.duke.edu 

Department of Mathematics, Duke University, Durham, NC 27708- 
0320, USA 



